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In thin magnetic films with strong perpendicular anisotropy and strong de- 
magnetizing field two ordered phases are possible. At low temperatures, 
perpendicularly oriented magnetic domains form a striped pattern. As tem- 
perature is increased the system can undergo a spin reorientation transition 
into a state with in-plane magnetization. Here we present Monte Carlo sim- 
ulations of such a magnetic film containing a periodic array of non-magnetic 
defects. We find that the presence of defects stabihzes parallel orienta- 
tion of stripes against thermal fluctuations at low temperatures. Above the 
spin reorientation temperature we find that defects favor perpendicular spin 
alignment and disrupt long range ordering of spin components parallel to 
the sample. This increases cone angle and reduces in plane correlations, 
leading to a reduction in the spontaneous magnetization. 

1 Introduction 

Quasi two dimensional ultra thin magnetic films engender a large area of 
theoretical and technical interest, due in part to the large variety of mag- 
netic properties that can be produced [T, "2] and their applications in data 



storage [3, 4J. For a sufficiently high ratio of dipole to exchange couphng 
strengths, the ground state of thin magnetic films can consist of magnetic 
stripe domains O El [7]. For ffims with a strong perpendicular anisotropy 
a second phase transition is possible, in which spins reorient resulting in a 
non zero magnetization parallel to the sample plane [51 [9l [TO]. 
There are a number of lithographic techniques that can be used to cre- 
ate nanometer scale magnetic structures [TTl[T2llT3l[n[[T5l[T6l[T7l[T8l[T9] . 
When compared with isotropic films, periodic magnetic nano structures have 
been shown to significantly alter macroscopic properties such as anisotropy 
[201 [21], magneto-resistance [20], cohesive field [22l |23] and spin reorienta- 
tion temperature [2il [25]. 

On the micro scale, magnetic stripe domains can appear with long range 
orientaional order [6l [26] or forming complex patterns [271 IB]- Nano 
scale patterning has been to used to create pinning sites for domain walls 
[281 [23 [301 1211 [321 [23]. When the period of pinning cites is comparable to 
the natural stripe width, long range orientational order can be stabilized 
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Theoretically these quasi two dimensional systems have been studied with 
a variety of methods. For two dimensional isotropic systems the prob- 
lem of melting is reasonably well understood [34j, in particular the spin 
reorientation transition and stripe melting have been studied analytically 
[351 [Ml [371 [38] and with computer simulation [391 SOI E]. Theoretically 
the problem of melting in two dimensional systems has been considered for 
the case of particles with a periodic potential [42j. The pinning of domain 
walls has been explored for both random [13] and periodic defects [^. Re- 
cently micro-magnetic computer simulations have explored the contribution 
of periodic defects and edge effects to magnetic reversal and hysteresis [45j . 
Here we perform Monte Carlo simulations on striped magnetic system in or- 
der understand the effect of periodic non magnetic defects on the thermally 
driven spin reorientation and stripe melting transitions. 
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2 Method 



The system is modeled as a two dimensional square array of Heisenberg 
spins E 5^, with lattice spacing a. 




where i and j represent two dimensional indexes, i — {ix^ iy), si — si- z and 
(. . . ) indicates the sum extends only over nearest neighbors. J, K and Cd 
represent the strength of the exchange coupling, perpendicular anisotropy 
and dipole coupling respectively. In order to introduce non-magnetic defects, 
some lattice sites are left empty. These defects are arranged as a regular 
square array with spacing Wd- The system is evaluated with Metropolis 
algorithm Monte Carlo. In order to approximate an infinite system, peri- 
odic boundary conditions are introduced. After the change of co-ordinates 
rnm = G + Pn'm' where = {px.Py) and p^.py G [0,L], dipole cou- 

pling is calculated over a series of replicas of the original system 06l ST] the 
dipole interaction and Monte Carlo steps are parallelized on a GPU using 
the stream processing method described in our previous paper [41j. Non 
magnetic sites si — are not updated. 

3 Results 

In order to create a periodic array of defects we select a system size of 
L — 64q: and defect spacing — 8a. The ratio /C = K/Cd was set as 
/C = 15 ensuring that the ground state was not canted (s^ • £ = ±1). The 
ratio of exchange to dipole coupling is selected to he J — J/Cd — 8.9 
giving a stripe width of Wg — 8a. At T = 0, when the system is ordered, we 
find that for the choice of parameters above, the lowest energy occurs when 
domain boundaries pass through magnetic defects (this minimizes the energy 
by replacing a high energy spin with a defect). The system is initiated in 
the ground state and Monte Carlo ensembles are generated disregarding the 
initial 10^ Monte Carlo steps to allow the system to equilibrate. A further 
5 X 10^ steps are taken with states recorded every 50 steps. Previously we 
determined that 50 steps allowed sufficient independence between ensemble 
configurations. In order to examine the effects of the defects results are 
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included from an identical simulation performed on a perfect lattice that we 
shall refer to as the isotropic case. 

When describing results we will refer to the normalized temperature T — 
ksTC^^. In Fig. [l] sample states are shown for low temperatures near 
to where orientational order is destroyed. In the absence of defects, as T 
is increased, the striped system initially undergoes roughening at the stripe 
boundaries. The roughened domain boundaries are associated with localized 
canting of the spins away from perpendicular alignment. As temperature is 
further increased the system undergoes bridging between stripes that leads 
to the destruction of long range orientational order. With the inclusion 
of defects the same general trends occur: stripe roughening followed by 
bridging and eventual destruction of long range order. However, the presence 
of defects stabilizes the striped order at higher temperatures. In addition, 
differences in morphology are observed. In the absence of defects the stripes 
display long wavelength undulations. In the presence of defects walls are 
pinned. Instead of long wavelength bending, fluctuations exist as roughening 
of the sections of wall between defects. Also, in contrast to the isotropic case, 
we observe that this initial roughening of stripes is not associated with the 
appearance of canted spins. 



Z X Y Z X Y 




Figure 1: Example of spin configurations near the loss of orientational order. 
Spins vales are indicated according the color scale shown below, defects are 
colored black. Columns from left to right: 5|, and followed by the 
same states in the presence of defects. Rows from top to bottom: T = 1, 
r = 2.5, r = 3, r = 3.25 and T = 3.5 
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In Fig. [2] the behavior of the two systems is shown at temperatures above 
the loss of orientational order. In both cases the system forms regions with 
spins canted towards in-plane alignment and the existence of long range 
order in the in-plane components. As temperature is increased the systems 
become increasing granular before reaching the paramagnetic limit. 
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Figure 2: Example of spin configurations at high temperature. Colors are 
used to indicate spins values and the loation of defects as in Fig. [l] Columns 
from left to right: 5|, and followed by the same states in the presence 
of defects. Rows from top to bottom: T = 3.75 , T = 4.5 ,T = 7 and T = 10 



3.1 Orientational Order Parameter 

In order to analyze the loss of orientational order we locate vertical and 
horizontal perpendicular domain walls by using nf^ and HQ 01] 

i,j h.n.n 

where v.n.n and h.n.n indicate that the sums should be taken over all pairs 
of spins which are nearest neighbors in the horizontal and vertical directions 
respectively. The orientational order is given by 

0. = {\ni-nl\/{ni + nl)) (3) 

With the inclusion of defects the sums in Eq. [2] are restricted to run over 
all pairs that are not defects. 
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In Fig. H] is plotted clS cl function of the normalized temperature T. 
At low T both systems display a striped array with smooth boundaries 




Figure 3: Orient ational Order Parameter as a function of T, black circles 
represent the uniform system, while gray dots represent the system in the 
presence of defects. The transition region is replotted with a finer T scale 
in the insert. 



corresponding to Oz = 1. We observe that, while the transition profile is 
similar, the presence of defects increases the transition temperature. 



3.2 In Plane Magnetic Order 

At high temperature when spins are no longer entirely perpendicular the 
system can display net magnetization parallel to the system plane. Letting 
= l/NJ2iSi and My = l/NJ2i4 (^ith TV = L^) , the in plane 
magnetization is 

M|| = ((M2 + M2)I). (4) 

In plane magnetic order can occur only when spins are canted away from 
the perpendicular alignment, in order to measure the degree of canting we 
use the cone angle 




where 9i is the zenith angle of the spin at site i. When calculating 77 and My 
in the presence of defects N is replaced with N' = N{w'^ — l)/w'^ to account 
for the fact that the defects don't contribute to the averages. In addition 
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to these two single site order parameters, we calculate correlation functions 
between different spins. Taking 9i and as the zenith and azimuthal angles 
of Si respectively, Si.sj = cos{9i) cos{9j) + sin(0^) sin{9j) cos{(j)i — (pj). Since 
we are interested in plane ordering we calculate 




bj)) for Si ^ and sj ^ 
otherwise 



(6) 



In order to calculate the correlation as a function of distance we define 



'ij (7) 



3 



Here 11 is the Heaviside Pi function (n(x) = 0(x+l/2)0(x— 1/2)) and A^^ is 

oiao oiao 
o[o\o oioi# 

Figure 4: The correlation between the circled sites on the left will not be 
equal to the correlation between the circled sites on the right due to the 
presence of a defect (indicated here by a black circle). 

the number of spins contained in the average Gi{r); = Il({rij —r) /a). 
Gi(r) is the average correlation of the spin at site i with spins at a radius r 
from i. In a spatially isotropic state one expects that Gi(r) should depend 
only on the separation between spins. Here the inclusion of periodic defects 
breaks the isotropy. The correlation between two spins separated by distance 
r will depend on the proximity of the spins to a defect. We define the 
following average; letting n = (L/wd) 

^^(^) = ^I1 Gi'ir) = E En((-i'i - r)/a)G,j 

(8) 

with 

= i + Wdax + Wdby for a, 6 G [1, n]. 

The meaning of this correlation function is elucidated in Fig. [5j Here Gi{r) 
calculates the correlation between a fixed spin at site i and spins at some 
fixed distance r. Since the system is not spatially isotropic we expect that 
Gi{r) will depend on i. Qi{r) averages the Gi{r) over all sites with equivalent 
proximity to their closest defect. In the absence of the symmetry breaking 
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Figure 5: Gi{r) calculates the average correlation between the spin at site 
i (red circle) with all spins within a fixed radius (upper left gray circle). 
Spins separated by integer combinations of the vectors w^x WdV will have 
equivalent proximity to their nearest defect. Gi(r) averages Gi(r) over these 
equivalent sites. 

defects and in a uniform phase Qi(r) is not dependent on i. 
In Fig. [6] the parallel magnetization is shown as a function of temperature. 
Here we see that in the presence of defects the magnetic ordering is sup- 
pressed, and that the peak magnetization is reduced by around 20%. In 
Fig. [7| we see that the degree of spin canting is reduced for 2 < T < 7, 
however this reduced spin canting is not sufficient to account for the re- 
duction in peak magnetization. In Fig. [sj Qi{r) is plotted for T = 4.5, 
corresponding to the peak in-plane magnetization. In the isotropic system 
Qi has slow monotonic decay with increasing distance between spins. For 
the non isotropic system Qi is calculated for two choices of i. The first choice 
is i as a nearest neighbor to a defect, in this case the correlation is strongly 
reduced for all r. The other choice is i at maximum distance from a defect, 
in this case the correlation is comparable to the isotropic case for small dis- 
tances. However the correlation strength decreases rapidly as r approaches 
r = Go: (the location of the closest defects). 

In addition to the reduction in correlation strength we observe a periodic 
structure in both the defect cases due to the periodic defect lattice. This 
effect is particularly strong for the case when the Gi is calculated for i a max- 
imum distance from defects, here certain values of r will correspond to the 
average including several defects simultaneously. In order to gain a measure 
of the average effect of defects on magnetic correlation we also simulated the 
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Figure 6: Parallel magnetization as a function of T, black circles represent 
the uniform system, while gray dots represent the system in the presence of 
defects. 



system with randomly located defects at T — 4.5. For this case we observe 
that Qi lies between the results for the ordered defects described above. 
In order to understand this reduced magnetic order close to the defects we 
show the spatial dependence of (77^) in Fig. 10 where we see that close to 
defects spins have a slightly increased average angle to the plane. In Fig. 
[9] this dependence is shown with randomly located defects. We observe the 
same local increase in cone angle near to defects. The effect is especially 
pronounced in the top left of the figure, where we observe a accumulation 
of defects associated with a region of significantly decreased canting. 
When Wd — ^ the average spacing between defects is large compared to the 
range of the local canting effect. In Fig. |8]we noted that, far from defects, 
short range correlations are comparable to those calculated for the isotropic 
system. In figures 12 and 13 we show the spacial dependence of (77^) with an 
increased defect density at T = 4.5. In the ordered case we have let w^ — ^ 
and we see that the cone angle is no longer correlated with defect location. 
In contrast, when the same number of defects are randomly spaced as in 
Fig. 13, clustering leaves areas where the cone angle remains small. 



In Fig. IT] we plot Qi{r) for the high density defects. Unlike the low density 
case Qi{r) does is not dependent on proximity to the ordered defects. We 
note also that the clustering effect means that the average short range cor- 
relation length is enhanced slightly when the defects are disordered. In all 
cases with high defect density the correlation length falls to zero at finite 
radius and so no in plane magnetization can form. In Table [l] we give the 
cone angle and magnetization for the cases described here and note that 
within the precision of the simulation, despite the differences in morphology 
and correlation length, the strength of the magnetic ordering is dependent 
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on the density rather than periodicity of the defects. 




Figure 7: 77 as a function of T, black circles represent the uniform system, 
while gray dots represent the system in the presence of defects. 
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Figure 8: Gi as a function of r (given here in units of a) at T = 4.5 for i; as a 
nearest neighbor to ordered defects (blue circles), i at a maximum distance 
from ordered defects (red squares), in the presence of random defects (yellow 
diamonds) and for a defect free case (green triangles). 



3.3 Fluctuations 

We now consider the effects of ordered low density defects on fluctuations 
as a function of T. We calculate the autocorrelation function 

aHX)^{{X-{X)f), (9) 

of the three order parameters My, and 77, which we denote ctm, and 



respectively. In Fig. 14 a^is plotted as a function of T, here we observe 
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Figure 9: (77^) in the presence of randomly spaced defects at T = 4.5, black 
squares indicate defects and the value of r]i is indicated on the scale below. 



that the peak fluctuations occur at a higher temperature in the presence of 
defects corresponding to the stabilization of the striped structure. 
ctm as a function of T is plotted in Fig. 
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In both cases the fluctuations 
display two peaks corresponding to the creation and destruction of in-plane 
order. The low temperature peak is shifted towards higher T in the presence 
of defects. 
In Fig. 
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shows the same trend for both the patterned and isotropic 
cases, a broad peak with maximum occuring at T = 4.5 corresponding to 



Table 1: Comparison of the properties of isotropic system, ordered defects, 
and random defects at T = 4.5. 



System 


M|| 


V 


Isotropic 


0.463 


0.529 


Low Density Defects 


0.363 


0.556 


Low Density Random Defects 


0.359 


0.556 


High Density Defects 


0.064 


0.622 


High Density Random Defects 


0.080 


0.614 
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Figure 10: (77^) in the presence of regularly spaced defects at T = 4.5, black 
squares indicate defects and the value of iji is indicated on the scale below. 
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Figure 11: as a function of r (given here in units of a) at T = 4.5 in the 
presence of a high density of defects for i; as a nearest neighbor to ordered 
defects (blue circles), i at a maximum distance from ordered defects (red 
squares) and in the presence of random defects (yellow diamonds). 

the peak in-plane magnetic order. 
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Figure 12: (77^) in the presence of closely spaced defects at T = 4.5, black 
squares indicate defects and the value of r]i is indicated on the scale below. 

4 Conclusions and Comments 

Monte Carlo simulations have been used to investigate the effects of non- 
magnetic defects on the stripe melting and spin reorientation transitions. We 
have shown that the inclusion of non magnetic defects with spacing compa- 
rable to the natural stripe width affects the melting of stripes by creating 
pinning sites for domain boundaries favoring parallel alignment of stripes. 
At higher temperatures the two measures of the spin reorientation transition 
(reduction of the cone angle and the appearance of in plane magnetization) 
are reduced. In particular there is a spatial dependence of cone angle and 
correlation strength on proximity to a defect. Recalling that dipole coupling 
favors in plane alignments of spins [39l [40], we surmise that the increased 
cone angle is due to the reduced dipole field near to the defects. The reduced 
dipole field increases the effective anisotropy near to the defects, suppressing 
canting away from perpendicular alignment. The increased 77 values reduce 
the size of the sm{0i) sm{0j) cos(0^ — 0j) term in s^.Sj, effectively reducing 
the exchange coupling of the x and y components of the spins, leading to 
the suppression of in-plane magnetization. 

Here we have restricted our attention to point defects. Recently Van de 
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Figure 13: (77^) in the presence of regularly spaced defects at T = 4.5, black 
squares indicate defects and the value of iji is indicated on the scale below. 

Wiele et al. have performed a temperature independent micro-magnetic 
simulation of magnetization reversal in a sample with square holes [45j. 
Here the defects have dimension comparable to the spacing between defects. 
They find that the local shape anisotropy of the holes significantly affects 
the reversal mechanism. In light of these calculations it would be inter- 
esting in the future to consider to consider the melting problem on larger 
lattices where the effects of changing the size and shape of defects could be 
investigated. 
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Figure 14: as a function of T for the isotropic case (I) and the patterned 
case (P). 
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